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This paper addresses the origin of the sihcate emission observed in PG QSOs, 
based on observations with the Spitzer Space Telescope. Scenarios based on 
the unified model suggest that sihcate emission in AGN arises mainly from the 
illuminated faces of the clouds in the torus at temperatures near sublimation. 
However, detections of silicate emission in Type 2 QSOs, and the estimated cool 
dust temperatures, argue for a more extended emission region. To investigate 
this issue we present the mid-infrared spectra of 23 QSOs. These spectra, and 
especially the silicate emission features at ~ 10 and ~ 18 fim, can be fitted using 
dusty narrow line region (NLR) models and a combination of black bodies. The 
bolometric luminosities of the QSOs allow us to derive the radial distances and 
covering factors for the silicate-emitting dust. The inferred radii are 100-200 
times larger than the dust sublimation radius, much larger than the expected 
dimensions of the inner torus. Our QSO mid-IR spectra are consistent with the 
bulk of the silicate dust emission arising from the dust in the innermost parts of 
the NLR. 

Subject headings: Infrared: galaxies - Galaxies: active - IR observations - 
Galaxies: silicate emission - quasars: unification model 



1. Introduction 



Unified schemes for active galactic nuclei (AGN) postulate an obscuring torus surround- 
ing an accreting super-massive black hole. Models predict that the infrared spectral energy 
distribution (SED) of the torus depends sensitively on its orientation, geometry and density 



distribution (e.g. Pier & Krolik 1992 


Granato & Danese 1994 


Efstathiou & Rowan-Robinson 


1995 


Granato et al. 


1997 


Nenkova et al. 


2002 


). In particular, the tori are predicted to ex- 



hibit prominent silicate dust features in either absorption or emission, depending on whether 
an AGN is viewed with the torus edge-on (Type 2) or face-on (Type 1). Previous failures to 
detect strong 9.7/im silicate emission in Type 1 AGN led to several proposed modifications 



of the unified model. For example modified grain size distributions have been assumed (Laor 



& Draine 1993 Maiolino et al. 2001) or a dumpiness of the torus invoked (Nenkova et al. 



& Draine 


1993 


2002 


)• 



The Spitzer Space Telescope (Spitzer), with its good mid-infrared (mid-IR) wavelength 



coverage and sensitivity, has drastically changed our view of this problem. Siebenmorgen 



et al. (2005) and Hao et al. (2005) reported the first Spitzer Infrared Spectrograph (IRS) 
detections of prominent silicate emission features in the mid-infrared spectra of several lu- 
minous quasars. Sturm et al. (2005) reported the first detection of 10 and 18 fim silicate 
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emission features in a low-luminosity LINER (NGC 3998). Comparison to the 10/18 fim 
feature ratio of optically thin emission from silicate dust at different temperatures suggests 
a modest temperature (~200K) of the emitting dust. 

The presence of prominent lO/im silicate emission features in AGN covering a broad 
range in luminosity could be taken as direct evidence for the existence of an obscuring torus. 
However, it is not at all clear whether this emission actually arises from the inner regions 
of a face-on torus. Depending on the size and composition of the grains, sublimation occurs 
between ~ 800 and 1500 K ( Kimura et al.|200"2 ) which is also the temperature range expected 
for silicate dust located near the hot inner torus wall. 

The lower temperature indicated for the emitting silicate dust can be interpreted as 
evidence for dust emission from regions located further away from the central heating source. 
Several arguments support such a scenario. Sihcate emission has also been detected in Type- 
2 QSOs ( Sturm et al. [2006 Teplitz et al.|2006 ), whereas for an edge-on view of the torus, one 
would expect to see silicate in absorption only. An extended emitting region, with dimensions 
much larger than the inner torus dimension, is fully consistent with this result. Broad-band 



10/xm imaging of several nearby AGN suggests extended mid-infrared continuum (Cameron 



eFaL]p93| [Tomono et al.|200H |Bock et al.|2'000{ [Radomski etaL][2003| [Packham et"aL][2005| . 



Efstathiou (2006) has modeled the sihcate emission of the Type-2 QSO IRASF10214-f4724 



(Teplitz et al. 2006) invoking extended NLR dust in addition to an AGN torus. Marshall et 



al. (2007) conclude that some of the optically thin warm emission in the QSO PG0804+761 



may emerge from regions beyond the torus and suggest clouds in the NLR as a possible origin 
of this emission. These arguments indicate that silicate emission may originate in extended 
regions (~ 100 pc dimension). 

The nature and location of the extended silicate-emitting region is not yet known. In this 
paper we explore in a quantitative way one plausible interpretation, namely the association 
of this cool dust with the NLR. To this end, we present fits of our QSO spectra with a 
superposition of NLR dust models and spectral components representing the innermost hot 
dust and the bulk of the inner structure emission (both related to the torus) as well as the 
large scale host emission. The fitted model and the bolometric luminosity of each source 
enable us to estimate the cool dust distance and its covering factor. We note that our models 
cannot exclude the possibility of a torus contribution to the observed silicate emission. 



In section 2A_ we describe our QSO sample. In section [3} we introduce the model com- 
ponents and detail our fitting procedure. We also describe how we estimate the silicate dust 



cloud distances and the related covering factors. In section [3^ we discuss the dependence on 
model parameters. In section |4| we present the results of the fits, which are then discussed 
in section |5} Finally section [6] summarizes our conclusions. 
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2.1. The PG QSO Sample 



The sources used in our study are part of the Spitzer spectroscopy component (PID 
3187, PI Veilleux) of the QSO/ULIRG evolutionary study (QUEST). The project and the 
sample are described in Schweitzer et al. ( 2006[ ) (hereafter paper 1). The QSO sample is 
largely drawn from that of Guyon (2002) and Guyon et al. ( |2006 ). It consists primarily of 
Palomar-Green (PG) QSOs ( [Schmidt fc Green|1983[ ) and covers the full ranges of bolometric 
luminosity ~ lO^^'^~^^L0(based on the absolute B band magnitude and the SED of Elvis 
et al. 1994), radio loudness, and infrared excess {h'L^{60fim) /LboI ~ 0.02-0.35) spanned by 



the local members of the PG QSO sample (see also Jester et al. 2005 for a recent view on 
selection effects in the PG sample). B2 2201+31A is not a PG QSO but is included in the 
sample because its B magnitude actually satisfies the PG QSO completeness criterion of 



Schmidt & Green (1983). The QUEST sample used in this paper includes 23 of 32 objects 



from the Guyon sample. We add one Palomar-Green object from the Guyon et al. sample 
previously observed by Spitzer (PG0050+124 = IZwl; Weedman et al. ( 2005[ )). Table [l] lists 
names and redshifts of all 23 QSOs in our sample, six of which are radio-loud. This sample 
covers a range from Mb = —21 to Mb = —26, with median Mb = —23.3. We assume a 
cosmology with IIo=70kms^^ Mpc^^, Qm = 0.3 and Q\ = 0.7 throughout the paper. 



2.2. Data reduction 

For the QSO sample, spectra were taken both at 5-14/xm in the low-resolution (SL short- 
low) mode and at 10-37/im in the high-resolution (SH short-high and LH long-high) modes 
of the IRS (Houck et al. 2004). Slit widths of 3.6" to 11.1" include much of the QSO hosts 
as well as the vicinity of the AGN. Our data reduction starts with the two-dimensional basic 
calibrated data (BCD) products provided by version 12 of the Spitzer pipeline reduction. 
We used our own IDL-based tools for removing outlying values for individual pixels and 



for sky subtraction, and SMART (Higdon et al. 2004) for extraction of the final spectra 



Small multiplicative corrections were applied to stitch together the individual orders of the 
low-resolution and high-resolution spectra, as well as additive corrections for residual offsets 
still found between the low-resolution spectra and the SH and LH high-resolution spectra 
after zodiacal light correction of the latter. Paper I discusses in greater detail the Spitzer 
IRS observations and our data reduction procedure. 
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Table 1. QSO sample 



Ubject 


z 


Ul 


log [L^ioo/ergs ) 


radio loud (L)/quiet (Q) 






Mpc 






(1) 


(2) 


(3) 


(4) 


(5) 


PG0026+129 


0.1420 


672 


44.66 





PG0050+124 (IZwl) 


0.0611 


274 


44.30 





PG0838+770 


0.1310 


615 


44.16 





PG0953+414 


0.2341 


1170 


45.11 





PGlOOl+054 


0.1605 


768 


44.25 





PG1004+130 


0.2400 


1203 


45.23 


L 


PG1116+215 


0.1765 


853 


45.13 





PG1126-041 (Mrkl298) 


0.0600 


269 


43.82 





PG1229+204 (Mrk771) 


0630 


283 


44.13 





PG1302-102 


0.2784 


1425 


45.17 


L 


PG1309+355 


0.1840 


893 


44.81 


L 


PG141 1+442 


0.0896 


410 


44.31 


Q 


PG1426+015 


0.0865 


395 


44.44 


Q 


PG1435-067 


0.1260 


590 


44.39 


Q 


PG1440+356 (Mrk478) 


0.0791 


359 


44.22 


Q 


PG1613+658 (Mrk876) 


0.1290 


605 


44.70 


Q 


PG1617+175 


0.1124 


522 


44.29 


Q 


PG1626+554 


0.1330 


626 


44.44 


Q 


PG1700+518 


0.2920 


1505 


45.68 


Q 


PG2214+139 (Mrk304) 


0.0658 


296 


44.40 


Q 


B2 2201+31A 


0.2950 


1553 


45.91 


L 


PG2251+113 


0.3255 


1706 


45.63 


L 


PG2349-014 


0.1740 


840 


45.21 


L 



Note. — 
Col. (1) — Source name. 
Col. (2) — Redshift. 

Col. (3) — Luminosity distance in Mpc for a Ho=70kms~"'^ Mpc"""^, ri,„ = 0.3 and f^A = 0.7 
cosmology 

Col. (4) — Continuum luminosity XLx at 5100A rest wavelength (from spectra by T. Boro- 
son, taken from [Netzer et ar| ( |2007[ )) 

Col. (5) — radio loudness for PC QSOs taken from Sanders (1989) and for B2 2201+31A 
from Hutchings & Neff (1992) 
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3. Modeling the PG QSO IRS Spectra 

3.1. Model Components 

We have developed a procedure to fit the spectra of our sources with different compo- 
nents that account for the presence of optically thick AGN-heated dust emission, silicate 
dust emission and reprocessed stellar emission from the dust in star-forming regions of the 
host galaxy. 

Following recent AGN models (e.g. Siebenmorgen fc Efstathiou||2005 Honig et al.||2006" 



Elitzur & Shlosman 2006 ), we assume a hot dust component located near the central accretion 
disc that is mainly responsible for the shortest wavelength infrared continuum emission. 
AGN models typically consider this dust to be located within a torus-like structure. The 



sublimation temperature for silicate dust grains ranges from 800 K to 1500 K (Kimura et 



al. 2002). This is also the temperature range expected for dust located at the inner surface 
of the torus with a typical sublimation radius of roughly Rgub — 0.5 ■ y/Lbom pc, where 



the bolometric luminosity, Li^om, is given in units of 10 erg/s (e.g., Barvainis 1987, 1992 



Granato et al. 


1997; 


Nenkova et al. 


2002 



th dependence on grain material and size). To 
account for this component, we introduce a hot black body spectrum with a temperature 
restricted to range between 1000 K and 1700 K. To take the temperature distribution within 
the torus into account as well, we introduce two additional black bodies as fit components 



(section 3.2) with temperatures between 150 K and 1000 K. In the PG QSOs, these three 
components represent the bulk of the torus emission. The black body temperatures are free 
parameters allowed to vary continuously between their specified limits. 



The fourth component represents cool dust. As explained in paper I and in [Netzer et 



al. (2007) (hereafter paper II), this component is assumed to be dominated by reprocessed 
stellar emission from star forming regions in the host galaxy. To account for this component 
we introduce a cool black body with temperature limited to 35-65 K, a typical range in 
starburst galaxies. In paper I we demonstrated that on average the 7.7PAH/FIR ratio in 
our PG QSOs is the same as in starburst dominated, local ULIRGs. This fact indicates 
similar properties for the star formation in both galaxy types. Hence the assumption of 
a cool black body (T=35-65 K) to account for the reprocessed stellar emission from star 
forming regions within the host galaxy is reasonable. We note that the IRS spectra do not 
include a substantial part of the FIR emission of the QSOs. This limitation tends to shift 
the inferred best-fit temperature of this component towards higher values. In addition to 
intrinsic variations of the temperature of starburst dust, this bias is the reason we introduce 
a limited temperature range for this component. 



To account for the contribution of star forming regions to the 4-40 /im spectra of our 
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sources, we use the ISO-SWS mid-IR spectrum of M82 (from Sturm et al. 2000) as a starburst 
template spectrum. In particular, this template is used to fit the broad polycyclic aromatic 
hydrocarbon (PAH) emission features found in the Spitzer-IKS range, which are clearly seen 
in some of our objects. These features are typical indicators for recent star formation. The 
ratio of PAH emission to FIR emission in galaxies is related to the average radiation field 



intensity, echoed in the 60/100 /xm flux ratio tracing the large grain temperature (Lu et al. 



2003). By scaling the contribution of the cool black body relative to the contribution of the 



M82 template, the fitting routine has the freedom to take such a difference between M82 and 
the PG QSOs into account. Using starburst-dominated ULIRG templates instead of M82 
would be an alternative possibility. To exclude ambiguities regarding the power source (AGN 
or starburst) of the template spectrum, we prefer here to use M82 which, is a well studied, 
typical starburst galaxy. We have also tested the fits using the starburst template of Brandl 
et al. (2006). We found no significant differences (in ip"^ or the contributing components) 
compared to the fits using the M82 template. This is due to the overall small contribution 
of the starburst template and the similarity between M82 and the Brandl template. The 
contribution of the M82 starburst template is found to be consistent with the PAH detections 
in paper I. 

An important new ingredient of the present work is dust emission from the NLR. To 
fit the silicate emission, which is visible in nearly all of our QSO spectra, we use dust NLR 



models based on those of Groves et al. (2006). In contrast to torus models, our NLR models 



account for cooler, optically thin dust that is located further away from the central AGN. The 
code used to generate these models, MAPPINGS ill, encompasses all dust-related processes, 
including stochastic heating, which allows individual small grains to reach high temperatures 
when heated by energetic photons. For the fits to the IRS spectra, we have removed the 
line emission in order to concentrate on the IR continuum and the silicate features. The 
main physics of the AGN IR modeling, including the dust composition and incident AGN 



spectrum, have been discussed in detail in Groves et al. (2006), and we briefly recapitulate 
the main parameters here. 



In Figure [T] we show the incident heating spectrum, which is a fit by Groves et al. (2006) 



of two power- laws with exponential cut offs to the observations of Elvis et al. (1994). We 



assume the gas abundances have their solar values ( Asplund et al.|[2005 ). Similarly, the dust 
depletions are based on local measurements ( Kimura et al.||2003 ) for these models (see Table 



1 in Groves et al. 2006). We note that the actual metallicities of the QSOs in our sample 



may be higher, and the actual dust depletions uncertain. The NLR models assume a mixture 



of silicaceous and graphitic dust (Laor & Draine 1993; Draine 1984 Weingartner & Draine 



2001 ) with a grain size distribution arising from a modified grain shattering profile, leading 
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to a smooth exponential cut-off in terms of tlie grain mass at botli ends of tlie distribution: 

dN{a)/da = fca-=^- \ ^ ^(,/._)3 , (1) 

witli k defined by tlie dust-to- gas ratio and a being the grain size. The minimum and 
maximum grain sizes are O.Ol/zm and 0.25/im respectively, for both grain types considered. 
PAHs are assumed to be destroyed within the harsh NLR environment. 

In contrast to previous constant-pressure cloud models, we have assumed a constant 
density structure of n(H) = lO^cm"^ for simplicity, and explored a range of 13 ionization 
parameter^ (incident fluxes) ranging from logf/ = 1 to —3 in steps of 0.3 dex. These values 
are typical for NLRs, with log U = 1.0 being an extreme case, and log U = —3 becoming too 
cool to contribute any significant silicate emission. The 13 ionization parameters correspond 
to the 13 NLR models (see Table |2]), that will be used later to fit the spectra. The density 
assumed is on the moderate to high end for NLRs, which are actually likely to be stratified 



in density (e.g.. Groves et al. 2004). However, the assumption of constant density allows us 
to define accurately the incident ionizing flux, the dominant parameter controlling the dust 
temperature, and hence the emission. If a lower density of n(H) = lO'^cm"^ is assumed, the 
general shape of the IR continuum appears the same for a given incident flux. However, the 
ionization parameter for a given incident flux increases by a factor of 10, resulting in stronger 
higher ionization lines. The assumed density has no influence on the estimated distances of 
clouds from the AGN and should be considered as an approximation. In section [5| we will 
discuss the NLR properties also taking line emission (e.g., [NeV]) into account. Using the 
line emission, we find that a density of n(H) = lO^cm"^ is an upper limit for the density of 
the line emitting region. 



As in the Groves et al. (2006) work, we have limited the column density of these models 
to A^(H) = 10^^-^cm~^. Fitting the QSO spectra with lower column densities (e.g., iV(H) = 
lO^^cm"^) leads to unreasonably large covering factors. Using a column density of iV(H) = 
]^g2i.5pj^-2 rggui^g in a median covering factor of ~0.16. Both the ratio of IR to optical/UV 
flux and Hubble Space Telescope [HST) imaging of AGN dKriss et al.|[l992| suggest a ~ 30% 



covering factor as an upper limit for the clouds in the NLR. Increasing the column depth 
of the NLR models will only affect our results minimally. By a column depth of N{H) = 
lO^^'^cm"^, approximately 70 % of the heating flux has already been absorbed and re-emitted 
by gas and dust, meaning that the preponderant shape of the IR emission is already in place. 
In addition, at larger column depths most of the dust is cool and emits most of its IR flux 



^The ionization parameter is a dimensionless ratio relating the ionizing flux to the particle density: 



U 
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Table 2: AGN IR Model Parameters 
Model Number \og{U) log( ^^cident Fhuc^ 



erg cm s 



1 


1.0 


5.56 


2 


0.6 


5.16 


3 


0.3 


4.86 


4 


0.0 


4.56 


5 


-0.3 


4.26 


6 


-0.6 


3.96 


7 


-1.0 


3.56 


8 


-1.3 


3.26 


9 


-1.6 


2.96 


10 


-2.0 


2.56 


11 


-2.3 


2.26 


12 


-2.6 


1.96 


13 


-3.0 


1.56 



Note. — The NLR models assume a constant density of 
n(H) = lO^cm^'^ for simplicity. 



12 1.2 0.12 0.012 0.0012 0.00012 

Qr 1 1 — ■ I I I I 1 1 ' — I — I I < I n| 1 — I — ■ I I I || ' — I — I I I ■■! 1 1 — I ■ I I I I 



10.0 IDO.O 

tiv (eV) 



Fig. 1. — Incident accretion disc spectrum from 



Groves et al. 



(2006). 
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Fig. 2. — Model NLR IR spectra from logf/ = 1.0 (top model) to —3.0 (bottom model) in 
steps of 0.3 dex. 
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at wavelengths greater than 30/im, outside the range of the IRS spectra. 

The resulting nebular IR spectra are presented in figure [2] They show the range in 
possible sihcate emission, as well as the relevant emission lines possible under these conditions 
assuming the above gas density. We note again that our fits did not include the NLR emission 
lines, but in section [5] we will discuss the implications for the NLR from a comparison of the 
modeled and measured line fiuxes. 



3.2. Fitting Procedure 

To fit the QSO spectra with the components presented in the previous section we devel- 
oped a template fitting tool based on the IDL routine MPFIT written by C.B. Markwardt]^ 
This routine uses the Levenberg-Marquardt technique and we use it here to minimize a 
modified x^-value (hereafter called ip"^) defined by 



,2 _ V^/ -^obs(A) --Fmod(A) -2 /r,x 



A 



where -Fobs (A) is the monochromatic QSO fiux and Fmod(A) is the model fiux made of the 
rebinned M82 template spectrum, the four black body spectra and the NLR dust model 
spectra. 



Thus 



mod 



where Bi are the Planck functions. All models and template spectra are binned to the 
IRS spectrum resolution for our sources. The first term in F^o^{X) sums over the four black 
bodies, while the second sum includes all possible templates (in our case M82 and one (or 
more) of the 13 NLR models). Each component of the model spectrum is normalized to 
its maximum within the IRS range and then scaled by a factor or hi. The normalization 
prevents large variations in these factors. The IRS spectra of our PG QSOs have been fitted 
using the IRS rest wavelength ranges plotted in Figure 5-8. These ranges vary from source 
to source depending on the intrinsic redshift of the object. 



^http: / / cow. physics .wise .edu / ~craigni / idl / idl.html 



Err{X) usually represents the 1-sigma uncertainty on Fobs(A), but in our case it is 
defined in a different way. From tests with different weighting schemes, we found that we 
best achieved a reasonable weighting and a good fit quality using: 



Err{X) 



Fpl{X) 



(4) 



The flux trend FpL is from a power law flt to the IRS source spectrum and docs not trace 
individual features. This choice has been made to achieve similar weights for the same 



sampling density in micrometers. To ensure that better sampled wavelength regions do not 
dominate the fit due merely to their finer sampling, we estimate AA using a neighborhood of 
five wavelength bins at each wavelength, and we make Err{\) proportional to \J~-^- The 5- 
37/xm IRS spectra have spectral resolution varying between ~60 at the short wavelength end 
and ~600 at the long wavelength end. Our weighting method ensures a good fit quality over 
the entire wavelength range while simultaneously being able to fit the PAH emission features 
located within lower fiux regions. To further increase the fit quality at shorter wavelength, 
Err{X) is proportional to a/A, which results in a stronger weighting at shorter wavelengths. 

As we focus only on the continuum emission, we have removed all emission lines from the 
original spectra and the NLR models. All lines were cut from the spectra, with the resulting 
gaps not considered by the fitting routine. Wc fit each source with one NLR model plus the 
four black bodies and the M82 starburst template and repeat this procedure for all NLR 
models. To obtain the best fitting model template (Fmod) for each of the NLR models, the 
contributions of the chosen NLR model and starburst template and the black bodies (oj) 
are allowed to vary, as well as the black body temperatures (Tj) within their specified limits 
(i.e., hot: 1700 K-1000 K, warmi and warmg: 1000 K-150 K and cool: 65 K- 35 K). Finally, 
for each of our PG QSOs, we determined the best fitting (minimum -0^) model template out 
of all the NLR models. 

We do not consider the black body spectra used in our fits to be physically meaningful. 
Instead, they are a reasonable, physically motivated approximation to a smooth underlying 
continuum, necessary to achieve (in combination with the NLR models) a good fit quality 
for our spectra. To account for the strong emission found at the shortest wavelengths within 
the IRS range, wc have to introduce a hot continuum component (the hottest black body) 
for all of our QSO spectra. The limited MIR wavelength range we cover prevents a robust 
determination of all parameters of this component, since the bulk of its emission is not 
covered by the IRS range. The same is true for the coolest black body. To achieve a good fit 
quality, we additionally introduce two black bodies with intermediate temperatures. Using 
just one black body of intermediate temperature leads to degeneracies within the lO/im 



relative flux deviations ( — 




) at different flux levels Fp^. A A is the local wavelength 
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region, with some sources having the ~ lOfim sihcate feature reproduced by the black body 
instead of the NLR model. To prevent this, we introduce a second intermediate black body, 
which allows for a smoother underlying continuum emission. Including even more black 
bodies was not necessary since a satisfying level of fit quality is reached with only the four 
black bodies. Since we do not interpret the black bodies as real physical components and 
have only a limited wavelength range, we will not discuss their properties. Note that for 
some sources, the fitting routine automatically minimizes the contributions of certain fit 
components. This effect can be most dramatically seen in PG 1626+554, shown in figure [5} 
for which only two components (one hot black body and one NLR model) contribute to the 
best fit model. 

Beyond our standard fitting procedure, we investigated several alternatives. First, we 
tested the fits using grey bodies proportional to A^" with a = 0, 1, 2 instead of black bodies. 
This change had no significant effect on the fit results {ip'^, contributing components) or their 
interpretation. Next we tested for a non-thermal contribution to the mid-IR spectra of our 
sources by extrapolating a power-law model (based on radio data from NED) to the MIR. We 
found that for the radio-quiet objects, the non-thermal contribution is a factor 1000-10000 
weaker than the MIR emission at 10 /xm and can be neglected. Three of the six radio-loud 
objects have an implied weak non-thermal contribution at 10 /xm at the level of 1-10% of the 
observed continuum. For the other three radio loud sources (PG1302-102, PG2251+113 and 
B2 2201+31 A), the non-thermal contribution could be significant. These sources have been 
fitted using a power law, with the two intermediate temperature black bodies removed. We 
find that including a power law non-thermal component does not change our result, since 
it provides fairly smooth underlying emission similar to that provided by our black bodies. 
The identification of our sources as radio loud/quiet is listed in table [l} 

The fitting procedure also allows us to use multiple NLR models simultaneously. A 



comparison between use of one and multiple NLR models is discussed in section 3.4.2 Ad- 



ditionally our fitting routine allows us to obscure single components by a foreground screen. 



The possible impact of extinction on our results is discussed in section 3.4.1 



3.3. Cloud Distances and Covering Factors 

To estimate the NLR cloud covering factors and distances (from the central source) we 
need to use the intrinsic AGN bolometric luminosity. In this context, we have to discriminate 
between primary radiation emitted by the accretion disk and secondary radiation like that 
reprocessed by the torus. For a detailed discussion of the bolometric correction, see paper II. 
We adopt here a bolometric correction defined by L^oi = 7 x L(5100), where L(5100) is the 
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5100A rest wavelength continuum luminosity (AL5100). As explained in paper II, this choice 
avoids a double counting problem, since it accounts only for primary radiation and does not 
include reprocessed dust emission. A larger bolometric correction would imply larger cloud 
distances and smaller covering factors. The values of L(5100) adopted here are based on 
the Boroson Sz Green (1992) observations, which used a relatively small aperture. Given the 
high luminosity AGN, the contributions of the host galaxies at this wavelength are negligible. 
Since the ground-based spectroscopy does not show any indication for a high extinction, the 
intrinsic L(5100) and the derived bolometric luminosities cannot be significant larger than 
the values used in this work. The values of L(5100) for the 23 sources are listed in table [l] 

Given the calculated integrated incident flux F^n (see table |2| for the best i/j"^ model, we 
can obtain the dust cloud distance, 



This in turn allows us to calculate the required NLR covering factor c given the NLR 
dust luminosity: 



(6) 

Rdust dust 

where F^lr^^^ is the NLR model flux estimated by integrating the fitted NLR model 
over the observed IRS wavelength range and -Fr^^^j is the model flux, integrated over the rest- 
frame wavelength range, at the distance Rdust from the central source. -Dl is the luminosity 
distance of the source (listed in table [T]). The distances and covering factors obtained in this 
way are listed in table [3j 



3.4. Fit components and uncertainties 

In this section, we discuss the impacts of extinction and the use of multiple NLR com- 
ponents on the fits. We also demonstrate the fit quality by fitting different NLR models to 
the spectrum of PG1004+130. 
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3.4- 1- The Effect of Extinction 



All of our Type 1 QSOs show blue continuum emission, indicating a low extinction 
[Ay < 2). Nevertheless, we tested the effect of a line of sight extinction on our fit results. 
For this we assumed a uniform screen with fixed extinction in the range Ay = — 5. This 
was applied to the hot and warm black bodies as well as to the NLR models using the 



extinction curve of Draine (2003). We then repeated the complete fitting procedure using 



the extinguished black body and model spectra. As expected, small {Ay < 2) extinction 
does not affect the fit results. Using Ay = 2.5 for some sources, the best fitting NLR model 
became slightly hotter. For the largest {Ay = 5) case, we find that, on average, the best 
NLR model becomes slightly hotter, resulting in a cloud distance lower by a factor of ~ 2. 
We assume this factor to be the upper limit on the uncertainty in the estimated distance due 
to extinction. All fits of the PG QSOs presented in the following sections assume Ay = 
to match the low observed extinction. For comparison, we also tested the fits with the 
extinction curve of Chiar & Tielens (2006). In general, the fit quality became worse (larger 



■j/'^-values) compared to fits with the curve of Draine (2003), while the best fitting NLR 
model did not change significantly. Due to the better fit quality, we use the extinction curve 



of Draine (2003) in the following. 



3.4-2. The number of NLR components 



The fitting procedure described in section 3^ assumes a single NLR component to 
infer the silicate dust distances and covering factors. Another possibility is contributions 
from several dusty clouds at different distances from the central source. To test how such a 
cloud distribution would influence our conclusions, we repeated the procedure using multiple 
NLR components. The fitting procedure then minimizes the contributions of those NLR 
components that do not increase the fit quality. The best fitting NLR models derived from 
this procedure have then been compared to the best fit single NLR model of the same source. 
We found that for most sources, the best fit single NLR model is also the main contributor 
in the multiple NLR model fit, but this is sometimes accompanied by a second model (with 
similar but a smaller contribution. Also, in some sources the cool black body is replaced 
by a cool NLR component whose temperature is too low to contribute to the silicate emission. 
In three sources, the hottest and strongest NLR contributions come from components slightly 
hotter than those preferred by the single NLR fits. Figure [3] shows a comparison between 
single and multiple NLR fits for PG0026+129. The fit quality as judged by the value of ip'^ 
does not change significantly (single model fit ip"^ = 0.0119 ; multiple model fit ip"^ = 0.0097). 
In summary, our conclusion of extended silicate emitting dust is independent of the number 



- 16 - 



of NLR components included in the fits. However, we cannot exclude the possibility of 
smaller contributions from other NLR components. 



3.4-3. Dependence on cloud distance 

In fig. |4]we demonstrate the quality of the fit for PG1004+130, which has been chosen as 
more or less representative for the whole sample. We compare the best-fit model with two fits 
using models that differ by a multiplicative factor of v^IO in Rdust- Inspection of the residuals 
shows that these two fits are significantly worse than the best fit model. The mismatch can 
most clearly be seen at ~10/im, where the neighboring models differ significantly from the 
data. In the lower right panel of fig. |4]we show ip'^ as a function of dust cloud distance for 
this QSO. It can be seen that ip"^ is a smooth function with a quite well defined minimum at 
a dust distance of Rdust=52 pc. 



4. Results 

In figure [5] we present the best fits for the complete QSO sample. In the top half of 
each diagram, we show the best-fit model with the observed spectrum. The data are shown 
in black and the model in red. Also shown are the individual model components: the NLR 
models (blue, solid curve), the M82 starburst template with PAHs (blue, dashed curve) and 
the four black bodies in green (dotted). The lower part of each figure shows the fit quality 
across the spectrum. As can be seen, the fit is satisfactory (within 20% of the flux density) 
for all sources. 

Table [3] lists the best fitting NLR models for each source in our sample and the estimated 
cloud distances and covering factors. The cloud distances range from 9 pc (PG 0050+124) 
to 263 pc (PG 2251+113), with a median value of 40 pc. This is larger than the expected 



dimension of the inner torus (a few pc; Jaffe et al. 2004). These large radii are compatible 



with an extended silicate emission region. Since we do not use any torus model in our fits, 
we are not able to exclude the torus as the origin of the silicate emission. However, we have 
demonstrated that the NLR is a viable alternative as a source or contributor of the silicate 
emission. Using the average value of Fin computed from our fits and EqnjS} the average 
scaling between radius and luminosity is given by: 

Rdust ^ 80{UomY^'pc (7) 



with the bolometric luminosity Lbowe given in units of 10^^ erg/s. In figure |6| we plot the 
distribution of the best fit models for the 23 sources. Each NLR model is related to a scaling 




Fig. 3. — Comparison of a single (left) and a multiple NLR model fit (right). In both fits 
similar cool NLR models produce the sihcate emission features. ( NLR model(s) (blue, solid 
curves); blackbodies (green, dotted curves); total model (red); observed spectrum (black) ) 
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PG 1004+130 PG 1004+130 




rest wavelength {fim) 

0.005 I I I I I 

10 100 1000 

Rd.=i in [pc] 

Fig. 4. — Comparison of fit results for PG 1004+130 using the best fit model {Rdust=52 pc) 
or neighboring models differing by a multiplicative factor of v^TO in radius (i?rf„st=16;164 
pc). Conversion from NLR component to radius is via Eq. 5. For this example we also 
present the corresponding black body temperatures; however we caution that their physical 
meaning is limited due to the limited wavelength range we cover. 
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relationship Rdust — x ■ (Lhoi^eY^'^'pc, where x is the scahng factor. The distribution peaks 
near x = 68, which corresponds to model # 5 with an incident flux of 10^''^^-^^ (table 
[2]). The model number (7^) is a running number (1-13) that identifles the NLR model. 
With increasing model number the NLR dust model becomes cooler and the scaling factor 
increases. In table [2| the model number is listed together with the ionization parameter and 
the incident flux assumed for the respective model. From flgure |6j it can also be seen that 
the hottest model {x = 15) with the largest incident flux, is never the best fit model. The 
same holds true for models cooler than model 7^ 9. This indicates that the range of incident 
fiuxes tried here is adequate to fit the dust properties in the NLRs of our 23 QSOs. The 
range in Rdust defined by the hottest and coolest models is indicated by the dashed lines in 
figure [7] 

The good fits we obtain support the idea that the narrow line region is a possible 
candidate for hosting the relatively cool silicate dust seen in emission. In figure [7} we plot 
the estimated dust cloud distances versus the bolometric luminosities of the PG QSOs. The 
solid line in this plot indicates the sublimation distance for the silicate grains calculated from 

~ 0.5(Lf,o;46)^''^ pc (section 3.1). Note that we assume here a direct exposure of the 



silicate grains by the central source. On average we find for the dust distance: 

Rdust ^ nORsub (8) 

The dashed lines in figure [7] indicate the distance range covered by our NLR dust models. 
The typical NLR sizes in our sources are about 5-10 times larger than the largest distances 
we estimate for the NLR clouds. We refer to section [5]3] for a discussion of the NLR size. 



We also compared our results to other known observations. The lower triangle in figure 
[7| is an upper limit on the torus radius estimated from MIR interferometric observations of 



the nucleus of the nearby Type 2 AGN, NGC 1068 ( |Jaffe et al.||2004D . The upper triangle 
indicates the expected silicate dust distance as estimated from our average scaling relation 
for this source. This example shows that the torus dimension is much smaller than what is 
expected for the NLR dust cloud distance in an AGN of this particular luminosity. We will 
discuss this further in section 5.1[ The crosses in figure [7] compare our estimate using the 



average scaling relation (upper cross) for the silicate dust cloud distance in the Type 1 QSO 



PG 0804+761 with an estimate of Marshall et al. (2007) (lower cross) using a decomposition 



technique similar to ours. For this comparison, we assumed the bolometric AGN luminosity 
given inlMarshall et al.|(|2007l) of Lhom = 0.48. For PG 0804+761, the two estimates of Rdust 



differ by a factor of 2.9 and are in good agreement considering the expected uncertainties 
and the general scatter of the cloud distances. 



The cloud covering factors are listed in table |3} They range from 0.09 (B2 2201+31A) 
to 0.50 (PG1309+355), excluding PGlOOl+054 which is the only source with a covering 
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factor larger than one (1.16). This implausibly large value could be the result of the source 
variation (e.g. S. Kaspi et al. (2000) , A. Wandel et al. (1999)) between the epochs of the 
optical and the Spitzer IRS observations. The median covering factor for the full sample, 
excluding PG 1001+054, is 0.16. In figure [8} we show the distribution of the covering factor 
excluding PGlOOl+054. 

Table |4] presents the contribution of the NLR dust components to the total spectrum 
in two ways. In the second column, we compare the relative contribution to the 6-25/im 
rest-frame wavelength range. In the third column, we list the NLR contribution to the 15/xm 
rest-wavelength continuum source flux. The table also shows the relative contribution of 
the M82 starburst template to the 6-25 /xm rest-frame wavelength range. The upper limits 
for this contribution have been calculated by scaling the M82 template to match the 3a 
upper limits for the 7.7 PAH feature (paper I). The contribution of the starburst template 
is strongly correlated with PAH strength. Most objects which show clear PAH features in 
their spectra have a strong M82 contribution (see figure [s] for an example). The median 
contribution of the NLR dust to the 6-25 /im restframe emission is 24 % and the full range 
is from 15 to 51 %. 

In the literature, the strength of the silicate emission in AGN spectra has usually been 
expressed as an equivalent width of the 9.7/im feature after interpolation of the underlying 
continuum emission. However, caution is required in using this method, since the possibility 
of silicate emission hidden by the presence of star formation may bias those estimates. For 
example, in the fits of PG1440+356 and PG1613+658 in figure [Sj it can be seen that this 
method would result in an underestimate of the silicate strength. The apparent smoothness 
of the spectrum around 10 /im is caused by the superposition of starformation related PAH 
features and a silicate dust contribution that is actually quite strong. To achieve a good 
estimate for the silicate strength in such composite sources, a decomposition procedure is 
desirable. 
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Fig. 5. — fit results: NLR model (blue, solid curve); M82 (blue, dashed curve); blackbodies 
(green, dotted curves); total model (red curve); observed spectrum (black curve) 
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5. Discussion 

We have successfully fitted the MIR Spitzer IRS spectra of 23 PG QSOs using cool NLR 
components to account for the silicate dust emission, in addition to components representing 
the underlying AGN continuum and the host galaxy's star formation. We have investigated 
the robustness of the fit results by showing that the fit quality has a clear optimum as a 
function of NLR model (i.e. distance from the central source, Fig. |4]), and by probing for the 
effect of additional foreground extinction even for these optically unobscured Type-1 QSOs. 
Fitting our spectra with multiple NLR components simultaneously to allow for a dust cloud 
distribution, we confirm the typical cloud distances derived using one single NLR component 
but note that we cannot exclude smaller contributions from different components. 

Our main result is that the silicate emission may arise in an extended region. The use of 
NLR models in conjunction with an underlying continuum is already suggestive of a model 
combining a central compact and hot dust continuum source with silicate-emitting clouds 
overlapping the NLR. Investigating in more detail the physical implications of the fit results 
can test this model as well as permit a discussion in the context of other scenarios of AGN 



dust emission, like more classical torus models or scenarios invoking disk winds (Koenigl & 



Kartje 1994 Elvis 2000 Elitzur & Shlosman 2006). 



5.1. Dust cloud distances 

The dust cloud distances we derive are on average 170 times larger than the dust sub- 
limation distance, which for our PG QSO sample reaches at maximum ~1 pc for the most 
luminous sources. From the distribution of the best fitting NLR models, we estimate an 
average scaling relation of Rdust — 80{Li,oi46Y^'^ pc for the cloud-source distance using the 
expression Lboi = 7 x L(5100). The resulting cloud-source distances range from ~10 pc to 
~260 pc. The estimated distances and the good fit quality using our NLR models suggest 
that the NLR may host the silicate dust seen in emission. Figure [6] shows that the best-fit 
model distribution covers quite a narrow model range. This narrow range is also the rea- 
son for the quite small scatter if we plot cloud distance versus bolometric AGN luminosity 
(Eqng, as seen in figure [7| and may be one of the most important findings. It refiects the 
fact that the range of models (and corresponding values of Fm) that fit the data best is a 
small sub-set of the total number of models. 



Our results are in good agreement with the estimates of Marshall et al. (2007) for the 
QSO PG 0804+761. These authors used a similar decomposition technique and estimated 
a silicate dust cloud distance of ~ 19pc, a factor ~ 2.9 below the ~55pc obtained from 
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Table 3. Results from NLR Model Fits 



Object 


Best Fit NLR Model 

-U.2 
# 


, (Incident Flux ) 
erg cm ^ s 


Cloud Distance 
pc 


Covering Factor^ 


PG0026+129 


6 


3.96 


54 


0.25 


PG0050+124 (IZwl) 


2 


5.16 


9 


0.42 


PG0838+770 


4 


4.56 


15 


0.27 


PG0953+414 


6 


3.96 


91 


0.23 


PGlOOl+054 


9 


2.96 


107 


1.16 


PG1004+130 


4 


4.56 


52 


0.33 


PG1116+215 


4 


4.56 


47 


0.14 


PG1126-041 (Mrkl298) 


6 


3.96 


21 


0.20 


PG1229+204 (Mrk771) 


5 


4.26 


21 


0.14 


PG1302-102 


6 


3.96 


98 


0.41 


PG1309+355 


6 


3.96 


64 


0.50 


PG1411+442 


4 


4.56 


18 


0.17 


PG1426+015 


5 


4.26 


30 


0.21 


PG1435-067 


3 


4.86 


14 


0.17 


PG1440+356 (Mrk478) 


4 


4.56 


16 


0.14 


PG1613+658 (Mrk876) 


5 


4.26 


40 


0.24 


PG1617+175 


3 


4.86 


13 


0.12 


PG1626+544 


2 


5.16 


11 


0.10 


PG1700+518 


5 


4.26 


124 


0.15 


PG2214+139 (Mrk304) 


4 


4.56 


20 


0.10 


B2 2201+31A 


6 


3.96 


229 


0.09 


PG2251+113 


7 


3.56 


263 


0.17 


PG2349-014 


8 


3.26 


229 


0.09 



Note. — 

(1) — estimated covering factor assuming Ho=70kms^^ Mpc~^, il,n = 0.3 and ilj^ — 0.7 

(2) - the model number is a running number (1-13) identifying the NLR model. Larger model numbers 
correspond to cooler dust (lower incident fluxes) . In table |2] the model numbers are listed together with the 
respective ionization parameters and the assumed incident fluxes. 
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Table 4. Contributions to the total flux 



Object 
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M82 6-25 lira 
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Pni99Q-l-904 rMrlrTTT^ 


94 
zt 


9^ 
zo 




PP 1 "^09 1 09 


90 
zy 


oO 


•Co 




44 


f^O 
OU 


<lo 


PG1411+442 


16 


19 


4 


PG1426+015 


23 


25 


5 


PG1435-067 


29 


35 


<8 


PG1440+356 (Mrk478) 


16 


15 


20 


PG1613+658 (Mrk876) 


20 


24 


13 


PG1617+175 


18 


17 


<5 


PG1626+544 


37 


38 


<9 


PG1700+518 


19 


24 


<4 


PG2214+139 (Mrk304) 


28 


34 


<3 


B2 2201+31A 


31 


47 


<3 


PG2251+113 


40 


54 


<5 


PG2349-014 


22 


26 


3 



Note. — 

— The total flux is measured across the same wavelength range as for the 
models. 

(1) — contribution of the NLR model to the total flux between 6 and 25 
^m (rest-frame) 

(2) — contribution of the NLR model to the total flux integrated between 
14.95 and 15.05 ;um (rest-frame) 

(3) — contribution of the M82 template to the total flux between 6 and 
25 /Ltm (rest-frame). Upper limits have been estimated from a conversion 
of the (3cr) 7.7PAH upper limits presented in paper L 
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model number 




1519 958 678 480 303 215 152 96 68 48 34 24 15 

scaling factor x 

Fig. 6. — Best fit NLR model distribution: Each model is related to a scaling law Rdust — 
X ■ (Lfco/46)^''^pc for the dust cloud distance, where x is the scaling factor and Li^om the 
bolometric AGN luminosity in units of 10^^ ergs~^. 
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10000.0 F 




Fig. 7. — NLR cloud distance versus L^oh compared to the dust sublimation distance Rgub 
(solid line). On average, Rdust is 170 times larger than Rgub- The dashed lines indicate 
the full distance range covered by the NLR models. The triangles indicate the upper limit 
for the torus size of NGC 1068 ( Jaffe et al.| 2004) (lower triangle) and the expected silicate 
dust distance at this AGN luminosity (upper triangle). The two crosses demonstrate the 
agreement for the silicate dust distance estimates for PC 0804+761 by Marshall et al. (2007) 
(lower cross) and based on our average scaling relation (upper cross) considering the general 
scatter of the estimated cloud distances. 
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Fig. 8. — Distribution of the derived covering factors (excluding PGlOOH-054, which has a 
derived covering factor of 1.16). Each bar has a width of 0.05 . 
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our average distance-luminosity scaling relation. As an example, in NGC1068, MIR VLT 
interferometry (Jaffe et al. 2004) finds compact mid- infrared emission components within 
r<2pc, which is much smaller than the expected distance of silicate emitting dust for an 
AGN of this luminosity (see figure [?]), and also much smaller than the additional larger scale 



emission seen in standard mid-infrared images of this AGN (Cameron et al. 1993; Bock et 



al. 2000 Tomono et al. 2001). This discrepancy offers additional support for the idea of 



a separation between the torus and the silicate emission region. Also, the tens of percent 
of extended mid-infrared emission suggested in this and other mid-infrared imaging studies 
of nearby Seyferts (Radomski et al. 2003 Packham et al. 2005) are in plausible agreement 
with the fraction of 6-25/im emission ascribed to the NLR by our fits to PG QSOs (with 
somewhat higher luminosity). We emphasize that our method does not assume any special 
geometry. Thus, while we are able to reach reasonable conclusions about the distance of the 
silicate dust, we cannot derive information about the angular distribution of the dust. 



5.2. Covering factors 

The median covering factor we derive for our QSO sample is 0.16 (excluding PGlOOl+054), 
with the overall distribution seen in figure [8] These values agree with other estimates based 
on narrow emission line imaging and equivalent line width measurements, which suggest NLR 



covering factors of < 30% (e.g. Kriss et al. 1992 Netzer & Laor 1993). We note that in our 



case, the derived covering factors depend on the assumed column density of our models. For 
this reason, the covering factor is uncertain due to the uncertainty of the real column density. 
Fitting our sources with lower column densities than that assumed here (iV(H) = lO^^'^cm"^) 
leads to too large covering factors. On the other hand, increasing the column will only result 
in a small reduction in covering factor, since most of the incoming flux is absorbed by the 
dust in a cloud of A^(H) = lO^^'^cm"^. The covering factor is basically determined by the 
emitted flux of the model. A larger column depth means higher dust absorption and stronger 
IR emission, which results in a lower covering factor. However, the relationship is not linear 
due to the change of the dust opacity with wavelength and the cooling of dust at greater 
depths. Source variability between the epochs of optical and our Spitzer observations may 
introduce another uncertainty. On average, this effect should cancel out, but it may explain 
the large covering factor (1.16) derived for PGlOOl+0054. Differences between single and 
multiple NLR models indicate further uncertainties in the covering factor, since the contri- 
bution of each NLR component is somewhat lower in the multiple NLR model fit (the flux 
in multiple NLR model fits can be shared by similar components). 

The distance to the center of the hottest dust is related to the AGN bolometric lumi- 
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nosity via our average scaling relation. If the cloud size is assumed to be constant, then its 
covering factor would decrease if the cloud is placed further away from the central source. 
Another possibility is that rather than having discrete clouds, the entire ionization cone 
could be filled with dust. In this case, the covering factor does not correlate with distance. 
To test if the covering factor is correlated with the NLR contribution to the total IR flux or 
the bolometric AGN luminosity, we tested for a correlation between the covering factor and 
the NLR contribution to the 15 fim continuum source flux as well as to the total AGN IR 
luminosity. No such correlation was found. The absence of a correlation may also reflect the 
uncertainty in our estimate of the covering factor 



5.3. Silicate emission and NLR properties 

The typical distance found here for the silicate emitting dust is some 170 times larger 
than the dust sublimation distance and hence much further than assumed in canonical torus 
models. However, it is also smaller than "typical" dimensions assumed for the NLR. For 
example, Bennert et al. ( |2002 ) have derived a luminosity dependent NLR dimension of 



roughly R(NLR)=2.1L([OIII])22^ kpc, where L([OIII])42 is the [OIII]5007A line luminosity 
in units of 10"*^ erg/s~^. For the sources in our sample, this translates to a typical NLR 
dimension of 2 kpc (the full range covered in our sample, using this relationship, is about 1-3 
kpc). Schmidt et al (2003) later investigated this relationship and found NLR dimensions 
that are about of factor 2 smaller for sources similar to the ones in our sample. These 



relationships were also discussed in Netzer et al. (2004) and shown to be inconsistent with 



dimensions derived for the most luminous AGN. Regardless of the exact value, it seems that 
the dimensions found here for the silicate emitting dust are 5-10 times smaller than the NLR 
size derived from the [OIII]5007A line luminosity. 

We have also looked at the line emission expected from gas clouds situated at the 
distance of the silicate emitting dust. Observed line fluxes will be published in a forthcoming 
paper (Veilleux et al. 2008, in preparation). Our continuum fitting procedure can only 
indicate the cloud distance. The line emission depends on the level of ionization of the gas 
(the ionization parameter) and hence on the gas density. While a detailed study of the NLR 
properties in QUEST QSOs is beyond the scope of the present work, we mention here several 
of the more important conclusions. 

• Assuming a density of about lO^cm"^, gives for the typical incident flux found here, 
an ionization parameter of about 10"°'^. For this ionization parameter, the high exci- 
tation lines of [NeV] at 14.3 and 24.3 fim can reach their observed luminosities, given 
the covering fractions we derive. However, the observed [NeV] 24/ [NeV] 14 line ratio 
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(typically found to be ~1 in our sample) suggests that this density is too high. 

• For a lower density gas, of about lO^ciW^, we get the required [NeV] line ratio but the 
implied ionization parameter is so high that the [NeVI]7.6 /xm line is predicted to be 
much stronger than observed. 

• In all cases of < lO^'^cm"^, the gas is too highly ionized to produce the strong 
observed lines of [OIV] 26 /im and [Nelll] 15.6 /im. These lines can be reproduced for 
much higher densities (lower ionization parameter) but, in this case, the [NeV] lines 
are predicted to be well below their observed luminosities. 

In short, there is no way to produce all the strong observed NLR lines from a single component 
situated at roughly the same distance as the silicate emitting dust. None of these conclusions 
are very sensitive to the gas metallicity, the exact shape of the ionizing continuum, the 
relative reddening of optical and IR emission lines, or the exact NLR geometry. 

Given this analysis we conclude that the region emitting the strong silicate features 
cannot represent the entire NLR. It can perhaps be associated with gas in the innermost 
NLR and the region emitting 'coronal' lines from highly ionized species. In this scenario the 
extended portion of the NLR that produces most of the [OIII]5007A emission, and probably 
also the intermediate ionization IR lines like [Nelll] 15.6 and [OIV]26 /im, is of much larger 
dimensions, perhaps by an order of magnitude or so. All this will be discussed in more detail 
in a forthcoming paper. 

While detailed models designed to fit the IR spectrum of the QUEST QSOs are not yet 



available, some examples of the above line relationships can be found in the recent Groves 



et al. (2006) study of the IR spectrum of dusty NLR clouds. 



5.4. Silicate emission from torus models 



Our finding of an extended silicate emitting region does not necessarily contradict torus 
models, since modified torus models (e.g. Nenkova et al. 2002 - Nenkova, M., Sirocky, M., 
Ivezic, Z., Elitzur, M., 2007, submitted to ApJ) allow for weak or absent silicate emission 
within a certain parameter range (e.g. a certain clump distribution). This is also true for 



some face-on torus models with smooth dust distributions (DuUemond & van Bemmel 2005 



Fritz et al. 2006 ) if certain parameters are assumed. A compact torus might be present and 



contribute strongly to the MIR continuum but less to the observed silicate emission. 



Many of the published torus models cited above predict some level of silicate emission 
for nearly face-on Type-1 configurations, but often with a higher ratio of the 10 and ISfim 
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silicate features than required by the apparently 'cool' dust found in our Spitzer spectra. In- 
depth comparison to the Spitzer spectra is needed and ongoing. To our knowledge, the first 
published result of this type is the fit of Spitzer-ohseived silicate emission features in two 



Type-1 AGN to a torus model presented by Fritz et al. (2006). These authors assume a flared 



disc with a continuous dust distribution, whose inner radius is defined by the sublimation 
radius of the dust, plus host galaxy emission (reprocessed stellar emission from dust in the 
host galaxy). They account for the different sublimation temperatures of silicate (Tsuf,=1000 
K) and graphite (T<j„f,=1500 K) dust grains which can result in a up to five times smaller 
minimum radius for the graphite grains. This difference results in a dust layer composed 
only of graphite grains, which attenuates the X-ray and UV emission from the central source 
and is taken into account when the authors compute the silicate minimum radius. 

Some of their fits (e.g., for PG1229+204 (Mrk771) and PG 2214+139 (Mrk304)) of 
Type-1 AGN, which are constrained only by broadband photometry, do not reproduce the 
silicate emission features observed with Spitzer (compare Figjsj). For PG 1229+204, their 
torus model shows silicate in absorption rather than in emission. The same is true for 
PG2214+139, which shows strong silicate emission features in its Spitzer spectrum. These 
discrepancies suggest that broad band data alone are not able to constrain properly the 
MIR dust properties, since the wavelength sampling is too low. In contrast, their fits to 
two Type-1 AGN with MIR Spitzer spectra (PG0804+761 and PGllOO+772 = 3C249.1) are 
satisfactory for both the sihcate features and the broad band SEDs. 



Two elements of the Fritz et al. (2006) models may contribute to this successful fit. First, 
the separate treatment of graphite and silicate sublimation moves the innermost silicate dust 
out in radius and behind some shielding graphite-only dust, resulting in cooler temperatures. 
Results will sensitively depend on the exact values assumed for the sublimation temperatures. 
Second, the two successful fits to Spitzer spectra are characterized by a very steep decrease 
of the total radial column density from the equatorial plane of the flared disc to its edge. 
While the column Nh is several lO^^cm"^ on the equator, it is about two orders of magnitude 
less near the surface. Both factors combined imply a noticeable covering factor by dust with 
a column of several lO^^cm"^ at radii of order 1.2-36pc, with the silicate grains setting in 
only at ~5pc. While this still falls short by a fair factor from the radii ~75pc implied by 
our scaling relation (for comparison we assume here the bolometric AGN luminosity used 
in Fritz et al. 2006), it suggests that part of the successful fit by this torus model may be 
due to a related concept: a moderate column of dust away from the equatorial plane of 
the torus and well outside the sublimation radius. Their fits to the MIR Spitzer spectra of 
PG0804+761 and especially PGllOO+772 are somewhat worse compared to our average fit 
quality, but we note that our assumed underlying continuum (black bodies) is less physical 
compared to a proper radiative transfer calculation. 
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A key challenge for torus models of silicate emission will be to explain silicate emission 



detected in some clear Type-2 AGNs (Sturm et al. 2006 Teplitz et al. 2006) 



The MIR wavelength range is very important for constraining model parameters, since 
the silicate emission features are very sensitive to dust temperatures, geometries, and density 
distributions. Salient AGN properties linked to the circum-AGN region include the Type 1/2 



ratio and the distribution of obscuring columns as determined from X-rays (e.g. Risaliti et 
ar||1999 ) or the optical/near-infrared obscuration of the BLR (e.g. Veilleux et al.||1997 Lutz 



et al. 2002). Models explaining silicate emission in terms of the torus alone need to fit those 



as well. A silicate emitting inner NLR of low optical depth, on the other hand, obviously has 
to be combined with obscuring structures that are plausibly related to the hot continuum 
in our decompositions. Finally, future high spatial resolution mid-infrared studies through 
8m-class diffraction limited imaging, and interferometry with improved sensitivity will be 
essential to remove remaining ambiguities between the contributions of torus and NLR to 
the silicate emission. 

Finally, while the match between our models and the source spectra is in general good, 
we find for some sources small deviations between the modeled and observed ~10 fim silicate 
peak position. Most clearly this can be seen in the fit of PG1626+554, where the observed 
silicate emission peaks somewhat redward from the modeled 9.7 fim silicate peak. This has 



already been reported earlier (e.g., Sturm et al. 2005) and may be explained by our physical 



assumptions, such as the grain size distribution and the exact chemistry of the silicate grains. 



6. Conclusions 

Scenarios based on the unified model suggest that silicate emission in AGN arises mainly 
from the illuminated faces of the clouds in the torus at temperatures near sublimation. How- 
ever, detections of silicate emission in Type 2 QSOs, and the estimated cool dust temper- 
atures, argue for an origin in a more extended region. To investigate this issue, we have 
presented the Spitzer-IRS spectra of 23 QSOs. We have matched physically-based models to 
the mid-infrared spectra and found that the silicate emission observed in these objects can 
be reproduced by emission from clouds, outside the central torus. This extended silicate- 
emitting region is possibly associated with the innermost NLR region or the intermediate 



dusty region proposed by Netzer & Laor (1993). 



The dust cloud distances found here scale with the AGN luminosity as Rdust ~ 80 ■ L^-^i^q 
pc, with the bolometric luminosity i^boZ46 given in units of 10^^ ergs~^. We have estimated 
the median distance of the dust cloud responsible for the silicate emission to be 40 pc, while 
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for individual sources distances up to 260 pc are possible. The smallest cloud distance is 9 
pc for PG 0050+124. The calculated covering factors for the dust clouds have a median of 
0.16 and are in agreement with an NLR origin of the silicate emission. Our models do not 
exclude the possibility of a torus contribution to the observed silicate emission, but rather 
emphasize the good agreement and perhaps the necessity of a larger-scale contribution to this 
emission. Finally, future high spatial resolution infrared observations and further crosschecks 
including the comparison to Type 1/2 ratios as well as distributions and individual values for 
the obscuring columns in X-rays and the near infrared are needed to resolve the remaining 
ambiguities between torus and NLR emission. 
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